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Abstract 


Shock-wave  effects  calculations  of  the 
Gasbuggy  gas  stimulation  experiment 
have  been  made  using  a  one-dimensional 
Lagrangian  code,  SOC.  The  detailed 
material  description  used  for  the  equation 
of  state  of  the  rocks  is  presented  along 


with  the  results  of  the  calculations.  Pre 
liminary  data  from  reentry  of  the  em¬ 
placement  hole  indicate  that  the  chimney 
height  and  available  time  of  arrival  data 
agree  with  predictions. 


Introduction 


The  purpose  of  this  paper  is  to  de¬ 
scribe  the  equation  of  state  used  and 
results  of  calculation  of  shock-wave 
effects  for  the  Gasbuggy  gas  stimulation 
experiment.  This  effort  was  motivated 
by  the  philosophy  that  the  successful 
industrial  application  of  nuclear  explo¬ 
sions  in  the  years  ahead  will  evolve  as 
the  result  of  a  systematic  scientific 
approach  to  understanding  each  experi¬ 
ment.  Such  a  philosophy  was  employed 
in  the  Gasbuggy  Event  as  a  precedent  for 
future  experiments. 

The  application  of  modern  computa¬ 
tional  techniques  to  underground  nuclear 
engineering  is  an  exciting,  yet  extremely 
challenging  new  field  of  research.  Both 
experimental  and  theoretical  methods 
are  employed  in  this  type  of  analysis. 

The  successful  application  of  these  tech¬ 
niques  has  opened  new  fields  of  scientifi¬ 
cally  directed  applications,  such  as  the 


use  of  nuclear  explosions  in  mining, 
excavation,  and  the  formation  of  under¬ 
ground  storage  cavities. 

The  growing  experience  in  this  field 
of  research  indicates  that  a  careful, 
representative  description  of  the  medium 
followed  by  calculations  (using  appro¬ 
priate  hydrodynamic  calculational  models) 
presents  a  much  improved  description  of 
shock  effects  when  compared  with  the 
often  used  scaling  techniques.1 

CALCULATIONAL  PROCEDURES 

The  models  currently  employed  for 

shock-wave  effects  calculations  are  SOC/ 

a  one-dimensional  Lagrangian  code  and 

TENSOR,  a  two-dimensional  Lagrangian 
2  3 

code.  *  These  codes  combine  the  equation 
of  motion  of  a  material  under  stress  with 
a  unique  description  of  that  material’s 
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behavior.  Figure  1  illustrates  the  cal- 
culational  process  for  establishing  a 
propagating  stress  field.  The  equation 
of  motion  is  obtained  by  combining  the 
conservation  of  mass,  linear  momentum, 
and  angular  momentum  equations  of 
continuum  mechanics.  This  equation,  in 
differenced  form,  provides  a  functional 


relation  between  stress  gradients  and 
accelerations.  References  2,  3  and  4 
discuss  the  differencing  methods  for  both 
the  equation  of  motion  and  the  strain 
tensor. 

In  this  calculational  process,  at  time  t, 
the  stress  is  known  in  each  mass 
element  in  the  medium.  The  stress 
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Fig.  1.  Feed  back  loop  for  stress  wave. 
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gradient  across  each  element  produces 
an  acceleration.  These  accelerations 
produce  new  velocities  which  in  turn  dis¬ 
place  (strain)  each  point  in  the  medium, 
and  cause  new  stresses.  Time  is  incre¬ 
mented  by  the  amount  used  to  produce 
displacement  from  acceleration.  This 
loop  is  repeated  until  the  response  of  the 
numerical  model  allows  a  prediction  of 
the  results  of  the  experiment. 

The  scheme  that  relates  strain  to 
stress  is  called  the  equation  of  state  of 


the  medium.  We  derive  this  description 
of  the  material  behavior  by  field  and 
laboratory  experiments  on  representative 
samples  of  the  shot  media.  Inherent  in 
the  equation  of  state  is  the  ability  to 
describe  elastic,  fractured,  plastic, 
liquid  and  gaseous  modes  of  behavior, 
and  also  to  allow  acceptable  transitions 
from  one  mode  to  another. 

In  order  to  obtain  a  meaningful  pre- 
shot  calculation  of  the  stress  induced 
effects  of  a  nuclear  explosion  on  the 
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surrounding  rock  media,  a  number  of 
logging  tests  and  laboratory  tests  on 
selected  core  samples  must  be  performed. 

The  logging  tests  pertinent  to  the 
equation  of  state  development  are: 

(1)  Density  log. 

(2)  Elastic  velocity  log. 

(a)  Compressional  velocity. 

(b)  Shear  velocity. 

Hopefully,  these  logs  will  permit  a 

judgment  of  both  the  layering  of  the 
medium  and  the  choice  of  core  for  labora¬ 
tory  testing. 

The  core  tests  required  are: 

(1)  Hydrostatic  compressibility  up  to 
40  kbar. 

(2)  Strength  tests. 

(a)  Triaxial  tests. 

(b)  Tension  tests. 

(3)  Hugoniot  elastic  limit. 

(4)  High-pressure  Hugoniot  data  and 
chemistry  for  those  rocks  near 
the  point  of  detonation. 

When  properly  interpreted,  these  tests 
give  the  strength,  compressibility,  and 
rigidity  modulus  of  a  material.  The 
source  description  is  obtained  by  calcu¬ 
lating  the  volume  of  rock  that  will  be 
shocked-vaporized  by  a  particular  yield. 
The  explosive  yield  is  then  distributed 
uniformly  within  the  vaporized  region  and 
the  stress  field  is  allowed  to  propagate 
from  this  source.  Reference  5  gives  the 
gas  equation  of  state  and  the  vaporized 
radius  for  a  number  of  rock  materials. 

THE  CHIMNEY  MODEL 

Strength  properties  of  granite  and 
dolomite  have  been  used  to  calculate  the 
extent  of  cracking  from  the  Hardhat  (5  kt. 


290  m  deep,  granite)  and  Handcar  (12  kt, 

402  m  deep,  dolomite)  Events. 

Figure  2  gives  the  strength  curves  for 
the  two  materials.  Figure  3  shows  the 
hydrostatic  compressibility  and  Poisson's 
ratio  for  each  material. 

As  mentioned  above,  the  strength, 
compressibility  and  Hugoniot  data  and 
Poisson's  ratio  (along  with  the  source 
description)  form  the  basic  input  for  the 
SOC  stress  wave  code.  The  results  of 
the  SOC  calculation  on  Hardhat  and  Hand- 
car  are  shown  in  Fig.  4.  It  is  interest¬ 
ing  that  the  calculated  extent  of  cracking 
for  both  Hardhat  and  Handcar  agrees  with 
the  observed  chimney  height  to  within 
14%.  The  calculated  and  measured  cavity 
radii  are  in  excellent  agreement  for  both 
experiments.  If  strength  of  material  and 
shock-induced  cracking  control  the 
chimney  formation,  then  it  is  no  longer 
surprising  that  a  12  kt  source  in  dolomite 
produces  a  smaller  chimney  than  a  5  kt 
source  in  granite. 

From  these  observations  it  is  proposed 
that  chimney  formation  should  extend 
above  the  cavity  a  distance  equal  to  or 
greater  than  the  shock-induced  cracking 
radius  predicted  by  the  SOC  code.  An 
exception  to  this  model  would  occur  when¬ 
ever  material  properties  were  such  that 
the  cavity  bulked  full  before  this  radius 
was  reached. 

GASBUGGY  EQUATION  OF  STATE 

The  important  formations  and  inter¬ 
face  positions  are  illustrated  in  Fig.  5 
for  the  Gasbuggy  experiment.  The  ex¬ 
periment  was  so  designed  that  a  complete 
equation  of  state  description  was  necessary 
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Fig.  4.  Calculations  for  Handcar  and  Hardhat. 
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Fig.  5.  Results  of  preliminary  Gasbuggy 
calculations. 
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for  the  Lewis  shale  formation.  The 
equation  of  state  of  the  Pictured  Cliffs 
sandstone  was  obtained  to  stress  levels 
of  200  kbar.  Only  limited  compression 
and  strength  data  was  required  for  the 
other  formations  since  they  were  exposed 
to  stress  levels  of  less  than  a  few  kilo- 
bars. 

Figure  6  shows  the  Hugoniot  data  that 
were  obtained  from  the  Lewis  shale  and 
the  Pictured  Cliffs  sandstone.  Figure  7 
compares  this  data  with  hydrostatic  com¬ 
pression  data  taken  by  D.  Stephens.  The 
offset  between  these  two  sets  of  data 
appears  (in  both  materials)  to  be 
consistent  with  the  triaxial  strength  data 
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shown  in  Fig.  8.  The  effective  strength 
as  measured  in  the  comparison  was  de¬ 
rived  from  the  expression 

°x  *  »m  +  4/3  K 

where  a  is  the  uniaxial  stress,  a  is 
x  m 

the  mean  stress,  and  K  is  the  strength  at 

a  given  mean  stress.  The  agreement  at 

the  two  different  strain  rates  indicates 

that  the  assumption  of  negligible  strain- 

rate  effects  in  the  use  of  statically 

determined  equation  of  state  data  is 

justified. 


The  equation  of  state  description  for 
the  Gasbuggy  materials  is  given  in 
Table  I.  These  data  have  been  chosen 
as  a  unique  and  representative  description 
of  the  media  involved  based  upon  the  tests 
mentioned  in  the  introduction.  Implied 
in  this  description  is  an  averaging 
process  inherent  in  the  logs  and  a  bias 
towards  the  more  competent  material 
used  in  the  laboratory  testing. 

CALCULATIONS 

The  input  equation  of  state  description 
presented  in  Table  I  was  used  in  two 
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Fig.  7.  The  hydrostats  and  calculated  Hugoniot. 

one -dimensional  calculations  assuming  predicted  a  cavity  radius  of  78  ft  and  a 

25  kt  of  energy  deposited  at  4250  ft.  The  crack  radius  extending  to  393  ft.  When 

first  calculation  ignored  all  interfaces  the  other  layers  were  included  in  the 

above  the  Pictured  Cliffs  formation  and  calculations,  the  Fruitland  coal  layer 
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Fig.  8.  Triaxial  strength  data  for  Lewis  shale  and  Pictured  Cliffs  sandstone. 


Table  I.  This  table  includes  the  parameters  that  are  used  in  conjunction  with  compres¬ 
sion  and  strength  data  to  generate  equation  of  state  descriptions  for  the  Gas- 

calculations.  The  heats  of  fusion  and  vaporization  that  are  given  are  for 

Si02. 


Material 

type 

Region 

thickness 

(ft) 

Lewis  shale 

Shot  point 
to  4202 

Pictured  Cliffs 
sandstone 

4202  to 
3915 

Fruitland  coal 

3915  to 
3882 

Fruitland 

shale 

3882  to 
3800 

Kirtland 

3800  to 

Comp.  Total  Compress- 

velocity  water  ibility 

Density  log  content  (mega- 

(g/cc)  (ft/sec)  (%  by  wt)  bars) 

2.61  13,700  3.5  0.191 


1.71  8,850 


2.49  11,700 


2.58  14.070 


Heat  of  Heat  of 

Poisson’s  vapor  fusion 

ratio  10l2erg/cc  10l2erg/cc 

0.15  4.759  X  10-1  1.259  X  10”1 


influenced  the  results  enough  to  stop  the 
cracking  at  that  boundary  334  ft  from  the 
detonation  point.  The  393  ft  cracking 
radius  (determined  in  the  first  calculation) 
is  a  measure  of  the  expected  lateral 
cracking  within  the  Pictured  Cliff  for¬ 
mation. 


TIME  —  msec 


Fig.  9.  Radial  stress  vs  time  for  the 
radius  at  10  meters. 


TIME  —  msec 


Fig.  11.  Radial  stress  vs  time  for  the 
radius  at  45  meters. 


Figures  9,  10,  11,  and  12  are  cal¬ 
culated  radial  stress  profiles  at  10,  20, 
45,  and  75  m  from  the  detonation 
point.  Figures  11  and  12  clearly  show 
the  elastic  precursor  propagating  ahead 
of  the  driving  stress-wave.  The  pre¬ 
cursor  stress  level  is  a  measure  of  the 


Fig.  10.  Radial  stress  vs  time  for  the 
radius  at  20  meters. 


T IME  —  msec 


Fig.  12.  Radial  stress  vs  tin^e  for  the 
radius  at  75  meters. 
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threshold  for  dynamic  failure  of  the  mate¬ 
rial.  Figure  9  shows  behavior  within  the 
Lewis  shale,  while  Figs.  10,  11,  and  12 
reveal  behavior  patterns  within  the 
Pictured  Cliffs  sandstone  formation. 

The  calculations  have  indicated  that, 
given  a  25  kt  source  at  4250  ft,  there  will 
be  material  failure  until  the  coal  seam  is 
reached  at  3915  ft.  This  highly  compres¬ 
sible  but  relatively  strong  material  will 
stop  the  cracks  from  propagating  above 
the  Pictured  Cliffs  formation.  As  the 
cavity  pressure  subsides  (and  collapse 
takes  place  to  fill  the  void)  a  competent 
coal  layer  should  prevent  collapse  above 
this  point.  However,  if  this  layer  proves 
unstable,  collapse  could  continue  until  a 
competent  formation  is  reached  (i.  e.,  into 
the  Fruitland  shale  formation). 


The  preliminary  data  that  are  now 
available  from  the  actual  experiment  in¬ 
dicates  that  the  energy  source  was  close 
to  the  expected  value  of  25  kt  (see 
Fig.  13).  Reentry  of  the  emplacement 
hole,  the  only  postshot  experiment  in 
progress,  has  provided  data  which  in¬ 
dicates  that  the  coal  layer  was  competent 
enough  to  halt  collapse  near  the  top  of 
the  Pictured  Cliffs  formation.  This  re¬ 
entry  drilling  encountered  a  6  ft  void  but 
no  radioactivity  at  a  depth  of  3858  ft. 
However,  at  a  depth  of  3907  ft,  a  9 -ft 
void  was  encountered  and  radioactive 
xenon  was  detected  at  the  surface  indicat¬ 
ing  that  the  chimney  had  been  located. 

These  preliminary  results  indicate 
that  an  accurate  equation  of  state  descrip¬ 
tion  of  an  experimental  environment  is 


Time  —msec 


Fig.  13.  Time  of  arrival  of  various  radial  distances. 


basic  to  reliable  predictive  calculations. 
Furthermore,  the  success  of  incorporat¬ 
ing  a  modern  failure  model  into  material 
behavior  has  dramatic  implications  when 


one  considers  the  broad  spectrum  of 
earth  material  strengths  that  are  en¬ 
countered  in  dealing  with  underground 
nuclear  applications. 
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